Scale invariance at the onset of turbulence in Couette flow 



o 

Oh- 
< 

0\ 



>v 

i 

q : 

o . 

C« . 
>> 

Oh' 



> 

in 

o 



s 



Liang Shi, x ' 2 >0 Marc Avila, 1 ' 3 ^ and Bjorn Hof 1 ' 4 '! 

1 Max Planck Institute for Dynamics and Self- Organization (MPIDS), 37077 Gottingen, Germany 

2 Institute of Geophysics, University of Gottingen, 37077 Gottingen, Germany 

3 Institute of Fluid Mechanics, Friedrich- Alexander- Universitat Erlangen-Niirnberg, 91058 Erlangen, Germany 

^Institute of Science and Technology Austria, 3400 Klosterneuburg, Austria 

(Dated: April 22, 2013) 

Laminar-turbulent intermittency is intrinsic to the transitional regime of a wide range of fluid 
flows including pipe, channel, boundary layer and Couette flow. In the latter turbulent spots can 
grow and form continuous stripes, yet in the stripe-normal direction they remain interspersed by 
laminar fluid. We carry out direct numerical simulations in a long narrow domain and observe 
that individual turbulent stripes are transient. In agreement with recent observations in pipe flow 
we find that turbulence becomes sustained at a distinct critical point once the spatial proliferation 
outweighs the inherent decaying process. By resolving the asymptotic size distributions close to 
criticality we can for the first time demonstrate scale invariance at the onset of turbulence. 



Turbulence often arises despite the linear stability of 
the laminar flow and the nature of this transition has 
remained unresolved for over a century pLj]. Prominent 
examples are pipe, channel and Couette flows and here 
finite amplitude perturbations can lead to sufficiently 
strong distortions of the base flow, that for large enough 
Reynolds numbers (Re = UL/is, characteristic length 
L, velocity U and kinematic viscosity v) causes transi- 
tion to turbulence. At moderate Re turbulent structures 
are localized and are commonly referred to as spots or 
stripes in Channel and Couette flow and puffs in pipe 
flow. Individual localized structures are of transient na- 
ture [2, 13J but they can also temporarily grow and seed 
new spots in their vicinity, overall resulting in complex 
spatio-temporal dynamics. As shown for pipe flow, this 
spreading rate increases with Re, eventually outweigh- 
ing the decay [4||. It has been argued that this mech- 
anism causes a phase transition in the thermodynamic 
limit from transient to sustained turbulence [5(. As Re 
surpasses the critical point, turbulence can overall sur- 
vive due to the spatial spreading. The observed dynam- 
ics bears resemblance to contact process such as directed 
percolation (DP) |fj|. In this analogy laminar flow corre- 
sponds to the absorbing passive state and turbulence to 
the active one. However, for pipe flow [4J the phase tran- 
sition has not been characterized directly but has only 
been inferred from the statistical behaviour (mean split- 
ting and decay times) of individual spots. The main dif- 
ficulty in pipe flow is that the relevant time scales of the 
relevant processes are extremely large, putting a study of 
correlation exponents presently beyond reach. Here we 
report a numerical (direct numerical simulation, DNS) 
study of another fundamental shear flow, plane Couette 
flow (pCf) where the fluid is sheared between two slid- 
ing plates (see Fig. [Ha)). As will be shown below, the 
timescales close to the critical point are much smaller 
than those in pipe flow. Hence size distributions and 
turbulent fractions can be determined close to criticality, 
which provides an unique opportunity to characterize the 



phase transition. 

Because of its unlimited width in the spanwise direc- 
tion (z' in Fig. [T]), pCf is spatially more complex than 
pipe flow. In order to simplify the spatial-temporal com- 
plexity and to enable us to resolve the long interaction 
times, we choose a slender computational domain in- 
clined with the flow streamwise direction, as shown in 
Fig. Hlb) ■ The idea of such a "tilted" computational do- 
main was originally introduced to reduce the computing 
cost by Barkley and Tuckerman pj. 



(a) Plane Couette Flow 




(b) 



(Le <<Lz) 




FIG. 1. Schematic of the computational domain for pCf. 
(a) normal box; (b) "tilted" slender box (white box in (a)). 
The stripe pattern in the background is a visualization of the 
streamwise vorticity field. 

Simulations are carried out using a numerical code 
for Taylor-Coucttc flow (flow between two rotating con- 
centric cylinders) J8J. The two cylinders are chosen to 
counter-rotate and a radius ratio (inner to outer cylinder 
radius) of -q — 0.993 has been selected so that curvature 



effects are negligible. As shown in an earlier study [9] 
Taylor- Couette flow approaches the pCf limit for these 
parameter settings and is linearly stable for all Re rel- 
evant to this study. By using the half gap distance h 
and th = h/U (here U is the magnitude of the bound- 
ary velocity.) as the length and time unit, the dimen- 
sionless imcompressible Navier-Stokes equations are nu- 
merically solved with a Petrov-Galcrkin pseudo-spectral 
method [8J]. We discretize the equations using Fourier 
modes in the in-plane directions and by modified Cheby- 
shev polynomials in the wall-normal direction. The 
Reynolds number is here defined as Re — Uh/v and the 
velocity field is decomposed into base flow and perturba- 
tion, u = Ub + u', where the perturbation part u' can 
be expressed as 



the probability that a stripe has not split up to a time t 

P(splitting at t^tf) =P- = l-(i-l)/N s ,i = 1, • • • ,N S . 

(2) 
As shown in Fig. EJa), the probability distributions of 
the splitting have exponential tails (excluding the initial 
formation period to)- This shows that stripe splitting is a 
memoryless process and we can therefore determine the 
mean time r s (Re) for a splitting to occur by the following 
exponential ansatz, 



P s (t)=cxp[-(t-t Q )/T s (Re)} 
with t s estimated by the sample mean 
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where ai nm are the spectral coefficients and v; rim are 
the Chebyshev polynomials. Periodic boundary condi- 
tions are imposed in the in-plane directions, while no- 
slip boundary conditions are applied in the wall-normal 
direction. The size of the computational domain is cho- 
sen sufficiently big in the z direction (L z x Lg x L r ) — 
(100hxl0hx2h) to study the evolution of isolated stripes 
whereas for the study of stripe interactions L z is in- 
creased to 960/i. The spatial resolution (number of grid 
points) for the smaller domain is (2L x 2N x (M + 1)) = 
(512 x 48 x 27), the adequacy of which is checked by 
the one-dimensional energy spectra and by the conver- 
gence study of the statistical lifetime and splitting time 
distributions. For the larger domain the resolution in z 
direction is increased to 6144. 

For the investigation of the decay and splitting of in- 
dividual localized turbulent stripes, uncorrelated veloc- 
ity fields of single stripes were used as initial conditions. 
These were generated by first simulating a fully turbulent 
flow at Re — 400, followed by a sudden reduction to the 
desired Reynolds number Re — 325. Since the decay and 
splitting are observed to depend strongly on the initial 
conditions, many realisations are necessary to determine 
the mean lifetimes and splitting times. The Reynolds 
numbers of interest are in the range of Re e [310,350]. 
At each Re, 300 realisations with different initial condi- 
tions are conducted and each simulation was run for a 
predefined time duration (cut-off time). If a decay or 
splitting occurred earlier, the run was terminated after 
this event. Let N s , N d and N c denote the number of 
splitting, decaying cases and the ones reaching the cut- 
off time, respectively. By sorting in increasing order all 
the final times accociated with spliting events (or split- 
ting time), we obtain a splitting time series {tf}fL 1 , with 



The sample mean is effectively the maximum likelihood 
estimator of the scale parameter t s [lOj . to has been de- 
termined the same way described in [3j. From Fig.[3Ja), 
we find that the mean splitting time (slope of the dis- 
tributions) decreases quickly with Re. The exact de- 
pendence of t s on Re is shown in Fig. [3] (dark square 
points). It turns out that the best fit for the simulation 
data is a superexponential function represented by the 
dark dashed line. Hence r s (Re) only approaches infinity 
asymptotically as Re decreases and consequently a non- 
zero splitting probability remains even for smaller Re. 

The same method is applied to obtain the probability 
distributions of the decay events (or lifetimes) , the result 
of which is plotted in Fig.[2][b). The lifetimes of localized 
turbulent stripes are also exponentially distributed, with 
the mean lifetime r scaling superexponentially with Re 
(red circles in Fig. [3J). Overall splitting and decay statis- 
tics show the same qualitative behaviour as pipe flow (see 
Ref. |J|), suggesting that the key physical processes are 
largely independent of the geometry and possibly apply 
to canonical shear flows. 

The intersection of the two curves t s ' (Re) in Fig. [3] 
fixes a distinct Reynolds number Re ~ 325 where the 
mean timescales of both processes are in balance, namely, 
r s (Re) = r d (Re). This value gives a lower bound for 
the critical Re c below which turbulence decays after suf- 
ficiently long time and it is also v ery close to earlier 
estimates of the critical point ll|, ll2| although those 



were carried out in domains with a much larger spanwise 
length. While characterisation of the spreading and de- 
cay processes for flows in a streamwise and spanwise long 
domains would require excessive computation time, the 
same line of argument would be applicable here. Again 
turbulence should become sustained once growth pro- 
cesses (streamwise and spanwise) outweigh the decay. 

In analogy to contact processes like DP the turbulent 
and laminar domains can be viewed respectively as the 
active and passive state. In order of turbulence to sur- 
vive in the system the splitting rate has to be larger than 
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FIG. 2. Probability distributions of single turbulent stripe 
at different Reynolds numbers, (a) Splitting times; (b) life- 
times, in the unit of Th = h/U. The dashed lines are the 
corresponding exponential curves from equation (J5J. Y-axis 
is in the logarithmic scale. 




FIG. 3. Scaling of the parameter t s ' with Re of the splitting 
(dark square) and decaying (red circle) of individual turbulent 
stripes. The error bars correspond to the 95% confidential 
intervals. The dashed lines are the superexponential curves 
to guide the eyes. 



the decay rate, otherwise turbulence dies out. The on- 
set of sustained turbulence may, in analogy to DP like 



systems, become sustained at a non-equilibrium phase 
transition |l3J . As a clear signature the spacing between 
active sites should become scale invariant close to the 
critical point and hence the passive regions have no char- 
acteristic length. 

To test this another set of simulations are performed 
and we extend the box size in the stripe-normal direc- 
tion e z to 960/i such that spatial correlations are taken 
into account. Through the simulations at three differ- 
ent Re (subcritical, critical and supercritical), the ob- 
served dynamics agree qualitatively with the DP model 
(see Fig. |4j: At Re — 300, the flow returns to the lam- 
inar state after sufficiently long time; At Re = 360, an 
initial single stripe quickly splits until it reaches a sta- 
tistically stationary state with a typical stripe spacing. 
At Re = 329.5 on the other hand turbulence persists 
but the stripe spacings change throughout, exploring all 
scales permitted in the given domain size. An example 
of stripe splitting is shown in Fig. 4(d). After an initial 
increase in width, the stripe breaks up into two segments 
of similar size that then continue to separate. It should 
be noted that this process differs from puff splitting in 
pipes where the new puff originates from a thin filament 
of vorticity that disconnects from the original puff (Fig. 
2 A in |4| ) . Furthermore in pipes splitting exclusively oc- 
curs in the downstream direction whereas in the present 
case no preferable splitting direction exists. 

For a quantitative evaluation of the laminar spacing, 
sizes of laminar gaps are measured by setting a cut-off 
value to the averaged vorticity, below which the flow is 
considered to be laminar, and the distributions are tested 
to be insensitive to the cut-off value (except for shifts of 
the absolute values). Data are sampled at the (quasi- 
stationary state, over a time interval of O(10 5 )Tf l - It 
is observed that the size distributions sufficiently above 
critical are exponential and hence possess a characteristic 
size. Close to the critical point the distributions follow 
a power law instead (Re — 329 in Fig. [5j, implying that 
there is no characteristic length. Size distributions hence 
indeed exhibit scale invariance, confirming that the inter- 
section point in Fig. [3] marks a phase transition and the 
onset of sustained turbulence. Furthermore the circum- 
stance that size distributions follow a power law close to 
criticality and become exponential at larger supercriti- 
cal values show that the resulting intermittent flows are 
intrinsically irregular and do not form fixed patterns as 
had been proposed previously 14|. 

In summary our investigation shows that the decay of 
turbulent stripes in Couette flow is a memoryless pro- 
cess and that individual stripes remain transient. Only 
through spatial proliferation turbulence can eventually 
become sustained. The splitting process that underlies 
this expansion of turbulence is memoryless. Overall this 
behaviour closely resembles the onset of turbulence in 
pipe flow and it is hence likely that these processes are 
generally responsible for the onset of turbulence in canon- 
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FIG. 4. Spatiotemporal diagrams at different regimes: (a) subcritical at Re = 300; (b) slightly above the critical point at 
Re — 329.5; (c) supercritical at Re — 360. (d) A close-up view of the splitting event at the moment indicated by the circle in 
(b). Snapshots are taken at the mid 6-z plane. Colormap shows the streamwise vorticity. 
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FIG. 5. Size distributions of laminar gaps at Re = 329 and 
Re = 340 in semi-logarithmic scale. The length in z direction 
is L z — 960/i. Inset is the plot in log-log scale. 



ical shear flows. In contrast to pipe flow here the relevant 
mean times of the decay and spreading processes are sev- 
eral orders of magnitude smaller which allows measure- 
ments close to criticality. The scale invariant distribu- 
tions observed show that the transition to turbulence is 
a non-equilibrium second order phase transition. While 
many aspects of this transition are analogues to directed 
percolation measurements of critical exponents and the 
universality class of the transition to turbulence would re- 
quire extremely long integration times and set challenges 
for future studies. 
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